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Abstract 

We investigate a stochastic model describing a column of grains in the jamming 
limit, in the presence of a low vibrational intensity. The key control parameter 
of the model, e, is a representation of granular shape, related to the reduced void 
space. Regularity and irregularity in grain shapes, respectively corresponding to ra- 
tional and irrational values of e, are shown to be centrally important in determining 
the statics and dynamics of the compaction process. 
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1 Introduction 



The study of slow dynamics in the jamming limit unifies the fields of granular com- 
paction [TJ |2] and glasses {3|. Key features of this involve frustration and hysteresis, 
among other complex phenomena [4], with the concomitant difficulty of modelling them 
in simple and physical ways. We present in the following a model of remarkable simplic- 
ity, which is nevertheless able to capture to a large extent the complex consequences of 
non-trivial interactions, even in one dimension. Issues that are probed include the effects 
of orientation, and thus shape, on packing in the jamming limit. 'Irregular' and 'regular' 
shapes of units (for example, grains) will be seen to have rather different consequences for 
compaction behaviour, when they are subjected to zero- and low-temperature dynamics. 

The present model is an extension, with interactions, of an earlier model of non- 
interacting grains, presented in For clarity, we summarise in what follows the com- 
monality and differences between the two. The previous model was two-dimensional. 
Each lattice site was occupied by an ordered grain (+), a disordered grain (— ), or a 
hole (0). It interpolated between a fluidised regime, where there were many holes on the 
lattice, and a jammed regime, where there were no holes anywhere on the lattice. The 
restricted dynamics in the jammed regime forbade migration of grains anywhere on the 
two-dimensional lattice, and in particular between columns. The result was expressible 
in terms of a column model of noninteracting grains, with a (trivial) ground state of 
completely ordered (+) grains. This, surprisingly, nevertheless exhibited some features 
usually associated with glassiness, such as slow dynamics and aging 

Clearly, disordered systems such as glasses or jammed granular media do not have 
ground states that are crystalline; equally, their attempts to reach their ground states 
are mediated by complex long-range interactions. The present model, already introduced 
in |Hj and investigated to some extent in jjj, represents an effort to make the jammed limit 
of the earlier model more realistic, by the inclusion of interactions. The column contains 
no holes. Each grain is either ordered (represented by a (+) Ising spin) or disordered 
(represented by a (— ) spin). However, and differently from the previous model, we take 
into account the effect of voids (holes that partially occupy a lattice site) which are 
associated with each disordered orientation of a grain. Thus, while each ordered grain 
fully occupies one unit of space, each disordered grain occupies e units of space, so that 
(1 — e) is a measure of the trapped void space. The net volume occupied by a disordered 
grain, or the corresponding void space, will depend on its shape; we see that e is thus 
a simple representation of granular shape. Also, and differently from the earlier model, 
the process of compaction is now no longer a simple relaxation into a completely ordered 
state: a given distribution of (+) and (— ) grains, as in nature, responds to externally 
imposed dynamics, by minimising void space locally, in the presence of disorder. The 
ground states so obtained resemble much more the random close-packed state found 
in granular systems |H] than the rather unrealistic crystalline state (all grains ordered) 
obtained before. 
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2 The model 




Our model is a fully directed model of interacting grains, where causality induces a 
directionality both in time and in space, as the orientation of a given grain only influences 
the grains below it, and at later times. Grains occupy all N sites of a column. As said 
above, we assume that they only assume two orientational states. We set o~ n — +1 
(resp. a n = — 1) if grain n is ordered (resp. disordered). A configuration of the system is 
uniquely defined by the orientation variables {cr n }. 

It is known that the response of jammed granular media to low- amplitude 
vibration is compacting; that is, grains rearrange themselves to locally maximise their 
packing fraction or, equivalently, minimise their void space. We model this by a local 
stepwise compacting dynamics: that is, a given grain orients itself to minimise void space 
locally, given the orientations of grains above itself. 

To be more specific, in the presence of a dimensionless vibration intensity T, we 
consider a stochastic dynamics, defined by the orientation-flipping rates 

= exp(-(A ri + /i n )/r), ( 
= exp(-(X n -h n )/T). 1 • 1 

In these expressions, h n and A n are, respectively, the local ordering field and the activation 
energy felt by grain n. These quantities are assumed to only depend on the orientations 
of grains above grain n. 

We make the further simplifying assumption that the activation energy X n does not 
depend on grain orientations at all. We write 

An = nl7£d y n, (2.2) 

where £d y n is defined to be the dynamical length [see (I2.6J) ]. Roughly speaking, £d yn is 
the depth of the nonequilibrium boundary layer: grains which are well within this length 
can order relatively freely in response to surface events, while grains much deeper relax 
only logarithmically in time 

The only dependence of the dynamics on orientations is via the ordering field h n , 
which determines the orientational response of grain n to orientations of grains above it. 
We choose to write the simple, linear formula 

hn = em~ - rra+ , (2.3) 

where m+ and m~ are the numbers of ordered and disordered grains above grain n: 

n—l i n—1 

m + = -£(l + ^), m- = -£(l-a fc ). (2.4) 
z k=i z k=i 

In spite of its simplicity, the dynamics defined by (|2.1|) - (|2.3|) can be shown to capture 
the compaction mechanism sketched in the Introduction, namely a local minimisation of 
the excess void space [TTj. In particular, a transition from the ordered to the disordered 
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state for grain n is hindered by the number of voids that are already above it. The whole 
picture will become clearer with the example of e — 1/2, discussed in section 3.2. 

In order to perform numerical Monte-Carlo simulations we will need a discrete-time 
formulation of the above rules. The flipping rates w n become flipping probabilities 




Pn 



l + exp(2/i n /r)' 



l + exp(-2/i n /r)' 



where the factor 



P n = exp(-A n /r) = exp(-72/£ dyn ) (2.6) 

describes the a priori exponential slowing down of the dynamics with depth n. 

Throughout the following, s, £d yn , and T (the latter being referred to as 'temperature') 
will be considered as three independent parameters of the model. In particular, e will 
not necessarily be restricted to the range < £ < 1, suggested by the interpretation of 
(1 — e) as the trapped void space. 

The zero-temperature statics and dynamics will crucially depend on whether e is ra- 
tional or irrational. Arguing that rational and irrational e corresponds to smooth/regular 
and rough/irregular shapes, respectively, this difference is to be expected, and will be 
discussed further. 



3 Zero-temperature statics 

As the dynamical rules (J2.1|) are fully directional, they clearly cannot obey detailed 
balance. The dynamics simplifies, however, in the r — ► limit j^j, where ()2.1|) yields 

w n { ►+) Joo if h n >0, 

—— = eX p{2K/T)^[ Q if ^ <0 (3.1) 

From a purely static viewpoint, ground states of the system can therefore be defined by 
the condition that the orientation of every grain is aligned along its local field, according 
to the deterministic equation: 

. J + if h n > 0, 

C7 n = srgn . f (3.2) 

provided h n ^ (see below). The condition ()3.2|) only involves the parameter e [see (|2.3j) ]. 
It is recursive, because of directionality, in that the right-hand side at depth n only in- 
volves the upper grains k = 1, . . . , n — 1. The uppermost orientation o~\ is left unspecified, 
as the corresponding local field vanishes identically. In the following, we assume for def- 
initeness that the uppermost grain is ordered: 

o v = +. (3.3) 
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It turns out that the zero-temperature rule ()3.2|) yields a rich ground-state structure, 
because of subtle commensurability and frustration effects. Our starting point is to 
observe that (|3.2|) implies 

( h n > =>• cr n = +, m + +1 = m + + l, m~ +1 = m~, /i n+i = /i n - 1, 

Un<0 cr„ = -, m + +1 = m+ , m~ +1 = m~ + 1, /i n+ i = h n + e. 

(3.4) 

The number and the nature of ground states depend on whether e is rational or irrational, 
which we consider separately below. 

3.1 Irrational e: unique quasiperiodic ground state 

Irrational values of e imply, in qualitative terms, a strong irregularity of grain shape. We 
refer to this as 'roughness'. Below, we demonstrate that the ground state of a packing of 
rough grains is unique, and optimally, but not maximally packed. We can visualise this 
as the interlocking of jutting edges to minimise, but not eliminate, voids. 

For irrational e, ()3.4|) implies recursively that all the local fields h n are non-zero, and 
that they lie in the bounded interval 

- 1 < h n < e. (3.5) 

Let us introduce the following superspace formalism. Consider the integers (m~,m+) 
as the co-ordinates of points on a square lattice. We thus obtain a broken, staircase- 
shaped line, starting as (mj",m+) = (0,0), (m^m^) = (0,1) [see (|3.3|) ]. etc. Vertical 
steps correspond to ordered (+) grains, whereas horizontal steps correspond to disordered 
(— ) grains. Equation ()3.5|) defines an oblique strip with slope e in the (m~,m + ) plane, 
which contains the entire broken line thus constructed (see Figure Q). 

A unique infinite configuration of grain orientations (i.e., a unique broken line) is 
thus generated. This configuration is quasiperiodic. Indeed the above construction is 
equivalent to the cut-and-project method of generating quasiperiodic tilings of the line, 
which has been extensively studied JTS] in the framework of quasicrystals. (Had we made 
the initial choice o~\ — — instead of ()3.3|) . we would have obtained the same quasiperiodic 
configuration, up to a permutation of the two uppermost grains.) We mention for further 
reference the following explicit expressions 1 for m+ and h n : 

m + = n-m- = l + Int((n-l)fi), K = -1 + Frac ^~ IM > (3.6) 



where the rotation number Q reads 



fi = e/(l + e). (3.7) 



1 Int(a;), the integer part of a real number x, is the largest integer less than or equal to x, and 
Frac(x) — x — Int(x) is the fractional part of x (0 < Frac(x) < 1). 
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An immediate consequence of (|3.6|) is that there are well-defined proportions of ordered 
and disordered grains in the ground state: 

f+ = n = e/{l + e), /_ = l-fi = l/(l + e). (3.8) 

This geometrical construction is illustrated in Figure Q for the most familiar irrational 
number, the inverse golden mean [T3*] : 

e = $ - 1 = 1/$, tt = 2 - $ = l/$ 2 , $ = (Vh + l)/2 « 1.618033. (3.9) 

The corresponding grain configuration is given by a Fibonacci sequence [T2"l I13j: 

Wn} = + + + - + + + - + + - + 




Figure 1: Geometrical construction of the quasiperiodic ground state of the model for the 
golden-mean slope (JEHU). The two ways of going around the first cell, marked with a circle, 
correspond to the two possible choices for the orientation of the uppermost grain. 

To recapitulate, at any point in the ground state, the ordering field h n will assume a 
finite, non-zero value, implying in its turn a given orientation of the grain that follows. 
The process is therefore fully deterministic, and the state obtained, unique. We will see 
that things will be rather different for rational values of e. 

3.2 Rational e: degenerate ground states 

Rational values of e imply, in qualitative terms, a regularity of grain shape or of void 
space, which we term 'smoothness'. We might here expect that regularly shaped grains 
could align themselves to fit into accumulated void space, in a given ground state con- 
figuration; this in fact happens, leading to states of perfect packing at various points of 



6 



the column. This in turn gives rise to the observed degeneracy of ground states to be 
discussed further below. 
For a rational e: 

e = p/q, tt = p/{p + q), (3.10) 

in irreducible form (p and q are mutual primes), some of the local fields h n generated 
by the recursion equations ()3.4|) vanish. The corresponding grain orientations a n remain 
unspecified. This means that grain n has a perfectly packed column above it, so that 
it is free to choose its orientation. For e = 1/2, for example, one can visualise that 
each disordered grain 'carries' a void half its size, so that units of perfect packing must 

be permutations of the triad H , where the two 'half voids from each of the (— ) 

grains are filled by the (+) grain. (Evidently this is a one-dimensional interpretation of 
packing, so that the serial existence of two half voids and a grain should be interpreted 
as the insertion of a grain into a full void in higher dimensions.) The dynamics, which is 

stepwise compacting, selects only two of these patterns, H and — I — (see the second 

line of Table 1). 

This feature of rational slopes is clearly visible on the geometrical construction. Fig- 
ure |2l corresponding to e = 2/3, shows that some of the lattice cells, marked with circles, 
are entirely contained in the closed strip (J3.5|) . Consider one such cell. The broken line 
enters the cell at its lower left corner and exits the cell at its upper right corner. It can 
go either counterclockwise, via the lower right corner, giving a n+ i = — , a n +2 = +, or 
clockwise, via the upper left corner, giving cr„ +1 = +, cr n+2 = — . Each marked cell thus 
generates a binary choice in the construction. This orientational indeterminacy occurs 
at points of perfect packing, such that n is a multiple of the period p + q, equal to the 
denominator of the rotation number Q. The model therefore has a non-zero ground-state 
entropy, or zero-temperature configurational entropy, £ = ln2/(p + q) per grain. Each 
ground state is a random sequence of two well-defined patterns of length p + q, each of 
them made of p ordered and q disordered ones, so that (|3.8|) still holds for each of the 
ground states. The patterns only differ by their first two orientations. The first cases 
are listed in Tabled Finally, the period p + q is formally infinite for an irrational slope. 
Accordingly, there is only one marked cell in Figure ^ for n = 0, corresponding to the 
fact that only the uppermost grain is unspecified. 



4 Zero-temperature dynamics 

Zero-temperature dynamics is a priori the canonical way of retrieving the ground states 
of a system. Here, the rule for zero-temperature dynamics is: 

cr„, ^ sign/i„, (4.1) 

according to (|3.1|) . with the definition (|2.3|) . We will find that irregular grains are able 
to retrieve their unique ground state, but that the degeneracy of the ground states for 
regularly shaped grains will make them impossible to retrieve. In the latter case, we find 
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Figure 2: Geometrical construction of the ground states of the model for the rational slope 
e = 2/3. The marked cells, entirely contained in the strip, are responsible for the non-zero 
configurational entropy. 



instead a steady state with non-trivial density fluctuations above the ground states, which 
recall the observed density fluctuations above the random close-packed state |H| ll()j. 

4.1 Irrational e, infinite £dyn : ballistic coarsening 

For irrational e, the rule ()4.1|) is always well-defined, as the local fields h n never vanish. 
We start with the situation where £d y n is infinite. We assume that the system is initially 
in a disordered state, where each grain is oriented at random: a n = ± with equal 
probabilities, except for the uppermost one, which is fixed according to (|3.3j) . 

The zero-temperature dynamics is observed to drive the system to its quasiperiodic 
ground state. This ordering propagates down the system from its top surface, via ballistic 
coarsening. At time t, the grain orientations have converged to their ground-state values, 
given by the above geometrical construction, in an upper layer whose depth is observed 
to grow linearly with time: 

L(t) ss Vt, (4.2) 

whereas the rest of the system is still nearly in its disordered initial state. 

This phenomenon is similar to phase ordering, as order propagates over a macroscopic 
length L(t) which grows forever. It is however different from usual coarsening, as the 
depth of the ordered region grows ballistically, vith a well-defined e-dependent ordering 
velocity V, instead of diffusively, or even more slowly [HJ. Figure H3 shows a plot of 
the inverse of the ordering velocity, measured in a numerical simulation, against e, for 
< e < 1. The ordering velocity obeys the symmetry property V{e) = V(l/e). It is 
observed to vary smoothly with e (although it is only defined for irrational e), and to 
diverge as V ~ \ je as e — > 0. 



8 



period p + q 


rot. number iz 


slope e 


P 


q 


pattern 1 


pattern 2 


2 


1/2 


1 


1 


l 


+ - 


- + 


3 


1/3 


1/2 


1 


2 


H 


- + - 


3 


2/3 


2 


2 


1 


H h 


- + + 


4 


1/4 


1/3 


1 


3 


+ 


- + 


4 


3/4 


3 


3 


1 


H h + 


- + + + 


5 


1/5 


1/4 


1 


4 


_|_ 


_ _|_ 


5 


2/5 


2/3 


2 


3 


_| 1_ _ 


1 h - 


5 


3/5 


3/2 


3 


2 


+ - + - + 


- + + - + 


5 


4/5 


4 


4 


1 


H h + + 


- + + + + 


6 


1/6 


1/5 


1 


5 


H 


1 


6 


5/6 


5 


5 


1 


+-++++ 


-+++++ 



Table 1: Patterns building up the random ground states for the first rational values of e. The 
second example with period 5 is illustrated in Figured 



4.2 Irrational e, finite £dyn- crossover to logarithmic coarsening 

For irrational e, in the situation where £dyn is finite, but large at the microscopic scale of 
a grain, the ballistic coarsening law ()4.2|) is to be modified as dL/dt ~ V exp(— L/^yn), 
taking the slowing down factor ()2.6)1 into account, hence 

L(t)«&ynln(l + W/£dyn). (4.3) 

Equation (|4.3|) exhibits a crossover between the ballistic law (|4.2j) for 1 Vt £<jyn, 
and the logarithmic coarsening law 

L{t) ^dynhlt, (4.4) 

already present in the model of non- interacting grains IE] • The dynamical length £d y n 
thus controls the spatial dependence of dynamical behaviour. In earlier work [5] it was 
shown to determine the extent to which order propagates down the column, in the glassy 
regime. This interpretation in terms of an ordered boundary layer continues to be valid 
in the present case: For an initially disordered state, the application of zero-temperature 
dynamics causes the quasiperiodic ground state to be recovered downwards from the free 
surface to a depth which grows ballistically with time. When L(t) becomes comparable 
with £d y n, the effects of the free surface begin to be damped, and in particular, for 
t 3> £,dyn/V, one recovers the logarithmic law ()4.4j) . widely associated with the slow 
dynamical relaxation of vibrated sand jU]. 

Equation (|4.3|) has been checked against the results of accurate numerical simulations, 
for the golden-mean slope. Figure 0] shows a scaling plot of numerical data for L(t) 
corresponding to £dyn = 50 and 100, together with the prediction ()4.3|) . with no adjustable 
parameter. The ordering velocity V ~ 2.58 is taken from the data of Figure El 
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Figure 3: Plot of the inverse ordering velocity 1/V of zero-temperature coarsening dynamics 
at infinite £dyn> against the irrational slope e, for < e < 1. 

4.3 Rational e, infinite £dyn- anomalous roughening 

We now turn to zero-temperature dynamics for rational e. The updating rule (|4.1j) is 
not always well-defined as it stands, as the local fields h n may now vanish. In such a 
circumstance, it is natural to choose the corresponding orientation at random: 

r + if h n > 0, 

cr„ , — > < ± with prob. 1/2 if h n = 0, (4.5) 

{- if h n <o. 

The zero-temperature dynamics defined in this way therefore keeps a stochastic compo- 
nent. We focus our attention onto the simplest rational case, i.e., e — 1. Equation (j2.3J) 
for the local fields reads 

n-1 
m=l 

We consider first the case where £d y n is infinite. We observe that the zero-temperature 
dynamics (J4.5|) does not drive the system to any of its degenerate dimerised ground 
states. The system rather shows a fast relaxation to a unique, non-trivial steady state, 
independent of the initial state. We now investigate this novel kind of zero-temperature 
steady state in some detail. 

Density fluctuations 

First of all, the local field h n has unbounded fluctuations in the steady state. Figure El 
shows that these fluctuations have a Gaussian distribution of width W n , at least deep 
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Figure 4: Scaling plot of L(f)/£d y n against i/£d y n for zero-temperature coarsening dynamics 
with the golden-mean slope. Symbols: numerical data. Full line: prediction (|4.3|) . with V = 
2.58. 



enough in the system (n ^> 1), except for a definite excess of small values of the local field: 
\h n \ ~ 1 <C W n . Figure |H1 (already shown in [7] and reproduced here for completeness) 
demonstrates that the local field variance grows as 

W* = {hl)^An*/\ (4.7) 

with A « 0.83. 

The exponent 2/3 of the anomalous roughening law (|4.7|) can be explained by means 
of the following local Markovian approximation. Assume that the local field h n obeys an 
effective Langevin equation of the form 

dh n /dt = -a n h n + r] n (t), (4.8) 

where r] n (t) is a white noise so that {i r lnif) r lnii')) = 2-D n 5(t — t'). Thus h n (t) is a Gaussian 
variable, whose steady-state width W n is given by the Einstein relation: 

W 2 n = (h 2 n ) = D n /a n . (4.9) 

The effective parameters a n and D n can be estimated as follows. For the deterministic 
part, ()4.5|) implies 

n-l 

d(h n )/dt =J2( a m- sign/i m ) « -(1 - Q n )(h n ), (4.10) 

m=l 
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Figure 5: Plot of the distribution of the local field h n for n as 1000. Histogram: numerical 
data (data for n = 999 and n = 1000 are mixed in order to avoid spurious parity effects). Full 
curve: Gaussian law with width VFiooo = 8.94. 



where the order parameter Q n is defined as 



Q n = (<7„sign/i n ). 



(4.11) 



The latter quantity will be shown below to fall off as n^ 1 ^ 3 [see (14. 22 J) ], implying a n « 1. 
The absence of divergence of the relaxation time r n = l/a n with n explains the observed 
fast relaxation to the steady state. As the fluctuating part is due to the second line 
of (|4.5jl . the strength of the noise D n reads, in some units, 



n-l 

D n nbY, Prob{/i m = 0} 

m=l 



n-l 

E 



2vr ^ W„ 



(4.12) 



assuming that the h n have a Gaussian distribution. Equation (|4.9J) yields 



W 2 w 

n 



n-l 

E 



i 



2vr ^ W„ 



(4.13) 



hence the power law (|4~7jl . with A = (96 2 /( 87r )) 1/3 - 

The anomalous roughening law (|4.7jl for the fluctuations of the ordering field h n is 
the most central feature of the zero-temperature steady state observed for rational e. 
It implies that, unless arranged by hand, the known ground states of a system of reg- 
ularly shaped objects (i.e., the crystalline ground states) will never be retrieved. On 
the contrary, the steady state will be one of density fluctuations above the ground state 
(which are related to fluctuations of the excess void space). The present model, to our 
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Figure 6: Log-log plot of W% = (/t 2 ) against depth n, for zero-temperature dynamics with 
e = 1. Full line: numerical data. Dashed line: fit to asymptotic behaviour, leading to (|4.7|) 
(after [Zj). 

knowledge, thus contains the first derivation of a possible source of density fluctuations 
in granular media [HIGH], which, here, arise quite naturally from the effect of shape. 

These anomalous fluctuations can be put in perspective with two different physical 
situations. First, as already underlined in [Zj, the power law (|4.7jl is reminiscent of the 
domain-growth mechanism in the low-temperature coarsening regime of the Ising chain 
with Kawasaki dynamics [15J. Second, using n to represent time in a random walk, 
the above results could be used to explain the ~ n 2 / 3 law recently observed in two- 
dimensional simulations of the effect of shape in particulate cages [HI] . 

Orientation and local field correlations 

If the grain orientations were statistically independent, i.e., uncorrelated, one would 
have the simple result (h^) = ne, while (14 .7|) implies that (/i 2 ) grows much more slowly 
than n. The orientational displacements of each grain are therefore fully anticorrelated. 
The anticorrelated orientational displacements are reminiscent of the bridge collapse seen 
in displacement-displacement correlations of strongly compacting grains JU]; grain ori- 
entational displacements in the direction of vibration were there seen to be strongly 
anticorrelated in jammed regions, as each grain tried to collapse into the void space 
trapped by its neighbours. Interestingly, correlations transverse to the shaking direction 
were ^U] found to be rather small, thus, in self-consistency terms justifying the choice 
of a column model in the present case. Once again, if we adopt a kinetic viewpoint and 
treat n as time in a random walk, these anticorrelations recall the temporal anticorrela- 
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tions observed in recent experiments investigating cage properties near the colloidal glass 
transition |17j . 

To be more specific, let us denote the orientation and local field correlation func- 
tions as 

Equation (|4.6|) implies 

m— 1 n— 1 

C m ,n ^ t ^ ] Ck,£i C-m^n Cm+l,n+l Cm+l,n Cm,n+1 "I - Cm,n; (4.15) 

fe=l 1=1 

and especially c n>n = C n+1>n+1 - 2C n>n+ i + C n>n = 1, so that C n , n - C n<n±1 « 1/2, and 
more generally 

C n ,„-C n ,n +fe ~ |fc|/2 (|fc|<n). (4.16) 
The power law (j4.7j) and the behaviour (|4.16j) can be combined into the scaling Ansatz 

(ft — rji \ 

where is a positive, even function, with a cusp at the origin of the form 

F{x) = 1 - \x\/2 + ■ ■ ■ (M<1). (4.18) 
As a consequence of (|4.15jl . the orientation correlations obey a similar scaling law: 

where -F(x) = d 2 jF/da; 2 is another positive, even function such that 

/+oo r+oo 
F(x)dx= xF(x)dx = l. (4.20) 
-oo JO 

The first of these sum rules confirms that spin fluctuations are asymptotically totally 
screened: J2 m ^n c n,n ~ _ c n,n = — 1- The scaling laws (|4.17|l and (|4.19|) are accurately 
confirmed by numerical data for C m n and c miU , whose scaling plots are respectively shown 
in Figures [7| and |H1 

A final consequence concerns the mixed correlation 

n— 1 

(cr n h n ) = c ™-™' ( 4 - 21 ) 

m=l 

for which the scaling results (|4.17j) . (|4.19|l yield (<J n h n ) ps 1/2. Scaling then implies that 
the order parameter defined in (|4.11j) falls off as Q n ~ l/W n , hence the estimate 

Q n ^an~ 1/3 . (4.22) 

This power-law decay is well confirmed by numerical data, shown in Figure El which yield 
a ps 0.44. 
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Figure 7: Scaling plot of the correlation function C m ,n of the local fields in the zero-temperature 
steady state with e = 1, demonstrating the validity of (|4.17|) . and showing a plot of the scaling 
function T . The full lines show the cusp behaviour (|4.18|) . 



All of the above adds up, for the case of regularly shaped grains, to a picture of grain 
orientations which are anticorrelated within a dynamical cluster jTO] whose size scales 
as n 2 / 3 ; grains outside such a cluster, are orient at ionally screened from each other, i.e., 
the screening length also goes as n 2//3 . (Correspondingly, from a kinetic viewpoint |16j . 
these results may be interpreted in terms of the time n 2 ! 3 spent by a walker bouncing 
back and forth between the walls of a cage, where his steps are consequently anticor- 
related one with the other.) Consistently, the order parameter Q n , proportional to the 
ratio of the screening length to the total length, goes as n 2 ^ 3 /n ~ n^ 1 ^ 3 . By contrast, 
when e is irrational, earlier orientations influence all successive ones, as the orientation 
correlations c mn do not decay to zero. The order parameter is Q n = 1 identically while, 
loosely speaking, the screening length scales as n. 

Entropy 

We now turn to the entropy of the steady state, denned by the usual Boltzmann formula 

S = -J2p(C)\np(C), (4.23) 
c 

where p(C) is the probability that the system is in the orientation configuration C in the 
steady state, and the sum runs over all the 2 n configurations C = {cr m } (m = 1, . . . , n) 
of a system of n grains. 
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Figure 8: Scaling plot of the orientation correlation function c mj „ for n / m in the zero- 
temperature steady state with e = 1, demonstrating the validity of ([4.19j) . and showing a plot 
of (minus) the scaling function F (after j7j). 



On the theoretical side, the entropy S can be estimated as follows, using the main 
feature of the zero-temperature steady state, i.e., the roughening law ([4.7)1 . Think of 
the depth n as a fictitious discrete time, and of the local field h n as the position of a 
random walker at time n. For a free lattice random walk of n steps, one has (/t 2 ) = n, 
and the entropy reads Sfl at = nln2, as all configurations are equally probable. Because 
{hn) = W% -C n, the entropy S of our random walk is reduced with respect to 5fl at . Let 

AS = S aat - S = nln2 - S (4.24) 

be the entropy reduction Consider first a strict constraint \h n \ < L. The probability 
that a random walk of n steps obeys this constraint is known to fall off exponentially, as 
V n ~ exp(— 7r 2 n/(2L 2 )). For a slowly time- dependent constraint \h n \ < L n , this estimate 
generalises to 

/ 7r 2 71 1 \ 
P„««p ~E7I ■ ( 4 - 25 ) 

\ Z m =l J^m ) 

With the assumption that the strict constraint \h n \ < W n and the weak constraint 
(/i 2 ) = W n generate similar entropy reductions for similar constraint profiles, we obtain 
the estimate 

n i 

AS = -lnP n ~ 5]_ 2 .~„V3. (4.26) 

m=l m 

We have evaluated the steady-state entropy S in a numerical simulation, using its 
definition (|4.23jl . by measuring the probabilities p(C) of all the configurations. As there 
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Figure 9: Plot of ra 1 / 3 times the order parameter Q n against n, for the zero-temperature steady 
state with e = 1, Symbols: numerical data. Full line: common limit value, yielding a ~ 0.44 
in m . 



are 2 n configurations for a system of n grains, the a priori statistical error only decays as 
(2 n /t) 1/2 . Reliable data are obtained in this way for t ~ 10 9 and n ~ 20. Figure ITUl shows 
a plot of the entropy reduction AS against n. The data show that AS* is small, at least 
for system sizes reachable by numerical simulations. For n = 12 (data of Figure ITTj) we 
have AS ~ 0.479. A reasonable semi- quantitative agreement with the estimate ()4.26|) is 
found: the fit shown in the plot suggests that (|4.26|) is affected by a logarithmic correction 
(which cannot be explained by the simple argument given above), with a small amplitude 
around 0.06. 

Figure ITT1 shows the normalised probabilities 2 12 p(C), plotted against the 2 12 = 4096 
configurations C of a column of 12 grains, sorted according to lexicographical order (read 
down the column). This plot exhibits a startlingly rugged structure on this microscopic 
scale: some configurations are clearly visited far more often than others. It turns out 
that the most visited configurations are the 2 6 = 64 ground states of the system (empty 
circles). We suggest that this behaviour is generic: i.e., the dynamics of compaction in the 
jammed state leads to a microscopic sampling of configuration space which is highly non- 
uniform, and reflects the structure of the ground states. In spite of this fine structure, the 
entropy reduction AS* ~ n 1 ^ 3 is subextensive, and therefore negligible with respect to the 
free entropy Sflat = n In 2, in qualitative agreement with Edwards' flatness hypothesis [TJ 
IT9~| . Our model thus provides a natural reconciliation between the intuitive perception 
that not all configurations can be equally visited in the jammed state where compaction 
is favoured, and the flatness hypothesis of Edwards, which is eminently sensible for 
macroscopically large systems. 
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with e = 1, defined in (|4.24|) . against n < 19. Symbols: numerical data. Full line: fit AS = 
(621nn + 53)10" 3 n 1 / 3 . 



4.4 Rational e, finite £dya : crossover to Brownian roughening 

In the case where £d y n is finite, the system still relaxes to a non-trivial steady state, which 
is qualitatively similar to that obtained for £d yn = oo, investigated above. 

At the quantitative level, the main effect of the fmiteness of £d y n is to induce a 
nontrivial profile of W%. In the regime where both n and £d y n are large, the following 
scaling law is observed 

C-Moc/W^n), (4.27) 

where (W%)oo is given by the anomalous roughening law (|4.7|) of the £d yn = oo steady 
state, which holds more generally for n <C £d yn , so that /(0) = 1. 

A qualitative understanding of the scaling function / can be obtained by generalising 
the above Markovian approximation. The expression (|4.12|) for the strength of the noise 
is readily replaced by 



V27T ~! W m 

For the deterministic part, (|4.5jl implies 

n-l 

d(h n )/dt =Y / (cr m - sign/i m ) e~ m ^. (4.29) 

m=l 

The right-hand side is not simply related to h n any more, so that a further level of 
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Figure 11: Plot of the normalised probabilities 2 12 p(C) of the configurations of a column of 
12 grains in the zero-temperature steady state with e = 1, against the configurations C in 
lexicographical order. The open circles mark the 2 6 = 64 ground-state configurations, which 
turn out to be the most probable (after [7]). 



approximation is needed. The most straightforward choice reads 

V£dyn ~ 



i n— 1 



n 



m=l 



n 



[1 - e ~ n/ ^y-). 



(4.30) 



Skipping the derivation, we mention that (|4.28|) . ()4.30|) imply ()4.27|h with 

,1/3 , , x \ 2/3 



X 



yf' 2 y" 1 ' 2 e-y 



dy 



i + x/i2 (x<i; 

K x 1 / 3 (a;>r 



(4.31) 



and K = 0.87732. In view of the crudeness of the above assumptions, (|4.31j) is only meant 
to provide a qualitative description of the scaling function /. Its asymptotic behaviour 
for x 3> 1 is, however, expected to yield the correct dependence 



AKi 



-1/3 
dyn 



n 



(4.32) 



of the width W n on n and £d yn - The profile of local fields is thus predicted to be Brownian 
for n 3> £d yn - Figure H2l shows a scaling plot of numerical data for the ratio W 2 /(W 2 )oo, 
against x = n/£d yn - A scaling law of the form 1)4.27)1 is clearly observed. The fitted curve 
is compatible with (j4.32j) . with K « 2.66. 
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Figure 12: Scaling plot of the ratio W%/(W%)oo against x = n/£d y n in the zero-temperature 
steady state with e = 1, illustrating the scaling law (|4.27j) . and showing a plot of the scaling 
function /. Symbols: numerical data. Curve: fit / = 1 + a((l + bx) 1 /^ — 1), with a = 5.63, 
b = 0.105, so that K = a b 1 / 3 = 2.66. 



5 Low-temperature dynamics 

We now turn to the investigation of the low-temperature dynamics of the model. Our 
main findings are the observation of intermittency in the position of the boundary layer; 
this has recently been observed in experiments of vibrated granular beds |2()j . 

We consider for simplicity the case of an infinite £d y n- If the slope e is irrational, the 
dynamical rule (j4.1|) is fully deterministic at zero temperature, so that a small non-zero 
temperature is expected to have drastic effects. To the contrary, for a rational slope e, 
the rule ()4.5|) is already stochastic at zero temperature, and indeed no interesting effect 
appears at a small non-zero temperature. 

We therefore focus our attention onto the case of an irrational slope e. We recall that 
the zero-temperature dynamics drives the system to its unique quasiperiodic ground 
state, where each orientation is aligned with its local field, according to ()3.2|) . For a low 
but non-zero temperature T, there will be mistakes, i.e., orientations o n = — sign/i n not 
aligned with their local field. Equation (J3.1|) suggests that the a priori probability of 
observing a mistake at site n scales as 

n(n) ss exp(-2|/i n |/r). (5.1) 

Hence the sites n such that the local field h n is relatively small in the ground state 
(\h n \ ~ r <C 1) will be nucleation sites for mistakes, and thus govern the low-temperature 
dynamics, in a sense that will become more precise. 
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The leading nucleation sites can be located as follows. Equation ()3.6|) shows that 
the local field h n is small when nVL is close to an integer m. The latter turns out to be 
m = m+. Indeed 

nVL = m + 5 h Tl = 5/(1 - O) (5.2) 

for 5 small enough (fl — 1 < S < O). The leading nucleation sites therefore correspond 
to the rational numbers m/n which are the closest to the irrational rotation number Q. 
This is a well-defined problem of Number Theory, referred to as the Diophantine approx- 
imation [T3] . 



5.1 The golden- mean slope 

Before we tackle the problem in general, we consider again for definiteness the golden- 
mean slope (|3.9|) . In this case, we are led to introduce the Fibonacci numbers F k [T^ ITB^. 
defined by the recursion formula 

F k = F k . 1 + F k . 2 (F = 0, F 1 = l). (5.3) 

We have alternatively 

F t = — LJ_. (5.4) 

The leading nucleation sites are the Fibonacci sites n = F k . We have m = m+ = F k ^2, 
m~ = F fc _i, and 

/ in = (-) fc $-( fc - 1 ), (5.5) 

so that 

n t =n W ~exp(-^). (5.6) 

We can therefore draw the following picture of low-temperature dynamics. Mistakes 
are nucleated at Fibonacci sites, according to a Poisson process. They are then advected 
with constant velocity V ~ 2.58, just as in the zero-temperature case. The system is 
ordered according to its quasiperiodic ground state in an upper layer (n < Af(t)), while 
the rest is disordered, somehow like the zero-temperature steady state for a rational slope. 
The depth Af(t) of the ordered layer, given by the position of the uppermost mistake, 
is a collective co-ordinate describing low-temperature dynamics. It evolves according 
to ballistic advection, i.e., J\f(ti) = Af(t ) + V(ti — t ), until it jumps backward to a 
smaller depth Af(ti) = F k , if another mistake is nucleated there. Figure HH1 shows a 
typical sawtooth plot of the instantaneous depth J\f(t), for a temperature T = 0.003. 
The leading nucleation sites are observed to be given by Fibonacci numbers. 

The system thus reaches a steady state, characterised by a finite ordering length (A/"), 
which diverges at low temperature, as mistakes become more and more rare. The law of 
this divergence can be predicted by the following argument. The most active nucleation 
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Figure 13: Plot of the instantaneous depth M(t) of the ordered layer, for the golden- mean slope 
at r = 0.003. Dashed lines: leading nucleation sites given by consecutive Fibonacci numbers 
(bottom to top: F n = 89, F 12 = 144, F 13 = 233) (after [7]). 



Fibonacci site is such that the nucleation time 1 /life is comparable to the advection time 
to the next nucleation site F fc+1 , i.e., (F k+ i — F k )/V ~ F k /(QV), hence the estimate 



1. 



(5.7) 



Indeed, less deep sites have too small nucleation rates, while the mistakes nucleated at 
deeper sites have little chance to be the uppermost ones. Equations (|5.6|) and (|5.7|) yield 



F k 

2$ §V 



(5. 



For T = 0.003, and for the Fibonacci sites shown in Figure the left-hand side of (|5.8J) 
respectively reads 0.56 for Fn = 89, 1.06 for F 12 = 144, and 1.94 for F 13 = 233. The 
estimate ()5.8|) therefore correctly predicts the observed fact that F i2 = 144 is the most 
active nucleation site at that temperature. 

The heuristic argument leading to ()5.8|) can be justified and made more precise by 
means of the results of Appendix A. The continuum approach used there is justified 
because the Fibonacci sites are more and more sparse. In the case of present interest, 
keeping only the Fibonacci sequence of leading nucleation sites, we obtain the predic- 
tion (jA.6|) for the ordering length, to be shown in Figure EJ 

For a low enough temperature T, the sum entering the right-hand side of f)A.6|) is 
sharply cutoff. It can indeed be argued that the term of order k in that sum is essentially 
F k _i for k < k*, while it is exponentially negligible for k > k* + 1, where k* = lnt(K), 
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and K is the real solution of ()5.8jh considered as a strict equality, with F K m $ K /\/5, 
according to (|5.4|) . We have therefore (A/ - ) ~ i 7 fc*+i, i-e., more explicitly, 

(AT)^F X ^. (5.9) 

The first factor of this expression, 

2$ / 1,2$ 



K 



V5T\\nT\ V l lnF l V5V\lnT 



ln-=- : + ••• , (5.10) 



shows that the ordering length obeys a linear divergence at low temperature, with explicit 
logarithmic corrections. The second factor, 

A K = ^-Fr-W (5.11) 

is a periodic function of its argument K « | In T | / ln$, with unit period, which oscillates 
between the bounds A max = $ and A m i n = 1. Oscillatory amplitudes are commonly 
observed in models related to self-similar structures 1TP : they originate in the discrete 
self-similarity of the underlying sequence. The oscillations of the asymptotic amplitude 
Ar, given in ([5. lip , are damped, except at extremely low temperature. Figure ITU shows 
a plot of numerical data for the product T(J\f), against |lnT|. These data are well 
described by the analytical prediction (|A.6|) . and lie within the bounds of the asymptotic 
estimate (|5.9j) - (j5.11j) . The oscillations become visible on the analytical curve for the 
lower temperatures (r < 10~ 4 ), which are not directly accessible to simulations. 

5.2 Other irrational slopes 

We now consider briefly the case of an arbitrary irrational slope e. The situation is rather 
similar to the phenomenon of hierarchical melting, observed at low temperature in some 
incommensurate modulated solids [2*2] . 

The nucleation sites can be determined as follows. The irrational rotation number Q 
can be written as an infinite continued-fraction expansion [T3] : 

0= — = [ai,a 2 ,a 3 ,...]. (5.12) 

ai H j 

a 2 + 



a 3 H 

The principal approximants of Q are the rationals 

tt k =Pk/qk, (5.13) 

whose numerators and denominators obey the same linear recursion 

Pk = a k Pk-i + Pk-2, Qk = a k qk~i + Qk-2, (5.14) 
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Figure 14: Plot of the product T(Af) against | In T | , for the golden-mean slope. Symbols: 
numerical data. Full line: analytical prediction (|A.6|) . Dashed lines: Extrema of the asymptotic 
result 1)5. 9 Jl . corresponding to A = -4 m ax (upper curve) and A = A m m (lower curve). 

with p = 1-\ = 0, pi = go = 1- The denominators build the leading sequence of 
nucleation sites. The rotation number Q also has secondary approximants 



for b = 1, . . . , ojt — 1 if cifc > 1. The denominators are the subleading nucleation sites. 

For the golden-mean slope, we have Q = 2 — $ = [2, 1, 1, 1, 1, . . .], so that the leading 
Fibonacci nucleation sites of section 15.11 are recovered, whereas there are no subleading 
nucleation sites. 

The most active nucleation site at low temperature can again be estimated by com- 
paring the nucleation time and the advection time. The ordering length (TV) is thus still 
predicted to diverge as 



at least for irrational numbers with typical Diophantine properties. Most irrational num- 
bers are typical in this respect. The presence of secondary approximants makes however 
the oscillation pattern of the amplitude ^4(lnT) more complex than a simple periodic 
function in general, in analogy with the low-temperature specific heat peaks induced by 
hierarchical melting [22] . 

In more physical terms, the ordering length (Af) defines the mean position of an 
intermittent boundary layer, separating an ordered state above it from a disordered 




(5.15) 




(5.16) 
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state below. This length is thus a kind of finite-temperature equivalent of the 'zero- 
temperature' length £ dyn . Both (Af) and £d y n retain the flavour of a boundary layer 
separating order from disorder. Within each of these boundary layers, the relaxation 
is fast, and based on single-particle relaxation, i.e., individual particles attaining their 
positions of optimal local packing ^UJ|2B| • The slow dynamics of cooperative relaxation 
only sets in for lengths beyond these, when the moves by which packing needs to be 
optimised become non-local. 

6 Discussion 

The work we have presented here concerns the effect of shape on compaction properties 
of grains. If irrational and rational values of e, the volume occupied by a disordered 
grain, are taken to correspond to irregular and regular grains, we see a distinct difference 
in behaviour between the two. 

While the uniqueness of packing of irregular grains in their ground state vis-a-vis 
the degeneracy of perfectly packed ground states for regularly shaped grains is in ac- 
cord with intuition, the effect of dynamics is more subtle. The fact that the perfectly 
packed and degenerate ground states of the regular grains are never retrieved by zero- 
temperature dynamics, leading instead to density fluctuations (that have been observed 
experimentally [9j) is rather subtle, as compared with the less perfect, unique and per- 
fectly retrievable ground state for the irregular grains. Clearly, a sharp distinction be- 
tween neighbouring rational and irrational values of e only makes sense for an infinitely 
deep system. For a finite column made of N grains, the distinction is rounded off by 
finite-size effects. In particular, the characteristic features of any 'large' rational e are no 
longer observed when the period p + q becomes larger than N. 

Zero-temperature steady-state density fluctuations for regular grains are subextensive: 
grain orientations are thus fully anticorrelated, reminiscent of dynamical heterogeneities 
in bridge collapse ^U] in strongly compacted granular media, as well as temporal anticor- 
relations in cages [THl HZ] . Also, while the macroscopic entropy of the steady state [TH] is 
approximately that of a fully disordered column, consistent with Edwards' 'flatness' hy- 
pothesis a microscopic examination of the configurations reveals a rugged landscape, 
with the most visited configurations corresponding to the ground states. As a matter 
of fact, the steady-state fluctuations of the local fields remain subextensive at any finite 
temperature. As a consequence, after a short transient regime, the steady-state propor- 
tions of ordered and disordered grains are solely determined by e, and given by (|3.8|) . for 
both rational and irrational values of e. 

Lastly, the low-temperature dynamics for irrational e leads to an intermittency of 
the surface layer. We would expect that for large enough temperatures, there would be 
little distinction between regular and irregular grains; our model however provides an 
interesting prediction of boundary-layer intermittency, which should be visible at low 
enough temperatures for irregularly shaped objects. 
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A Distribution of the depth M of the upper layer 



This Appendix is devoted to the depth N(t) of the ordered layer for low-temperature 
dynamics in the irrational case. Our main goal is to derive the stationary distribution 
of A/", for the effective dynamics described in section 15.11 

For convenience we use a continuous formalism, treating M as a real variable. Let 
it(x) dx be the given nucleation rate per unit time between x and x + dx, and p(x, t) dx 
be the probability of finding the depth M between x and x + dx at time t. The unknown 
probability distribution function p(x, t) obeys the rate equation 

\L + V ~IL) p{x ' t} = 71 {x)P{x ' f) ~ p{x ' t)u{x) = £ {n(x)P& *)) , (A.i) 

with the notations 

P(x,t) = £°p(y,t)dy, p(x,t) = - 9P ^ t \ U(x) = jT ir(y) dy, tt(x) = ^M. 

(A.2) 

Indeed, the left-hand side of (lA.lj) is the usual covariant derivative, whose convective term 
involves the drift velocity V. The middle side represents the evolution due to nucleation 
events, with the first (gain) term originating in nucleation at depth x, and the second 
(loss) term originating in nucleation at depth y < x. 

The stationary (time-independent) solution Pstat(^) of (|A.1|) is such that 

V Pstat (x) = -V dPst ; t{x) = n(x)P stat (x). (A.3) 
dx 

This separable differential equation easily yields the results 

P stat (x) = exp (-1 £ U(y) dy^j , p stat (x) = 5p exp (-± £ U(y) dy^j , (A.4) 
and especially 

(A/) = J exp i^-y J q n(y)dy^ dx = J exp (-— (x - y)ir(y) dy) dx. (A.5) 

These expressions hold for an arbitrary distribution of nucleation rates. 

In the case of interest in section 15.11 taking into account the leading sequence of 
Fibonacci sites F k , with nucleation rates H^, we obtain 



OO T T 

(AO = £ T- {^ Bk/V - e~ Bk+l/V ) > ( A -6) 

k=0 Ak 



with 

A k = £ n,, B k = £(F fe - F e )U e . (A.7) 
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